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We show that the addition spectra of semiconductor quantum dots in the presence of magnetic 
field can be studied through a theoretical scheme that allows an accurate and practical treatment 
of the single particle states and electron-electron interaction up to large numbers of electrons. The 
calculated addition spectra exhibit the typical structures of Hund-like shell filling, and account for 
recent experimental findings. A full three dimensional description of Coulomb interaction is found to 
be essential for predicting the conductance characteristics of few-electron semiconductor structures. 



In semiconductor quantum dots (QDs) Coulomb correlation effects are expected to become more and more pro- 
nounced as their spatial confinement is increased. This has been sofar one of the reasons of interest for such systems: 
they display the rich shell structure typical of atomic physics, but in addition they offer the possibihty to tailor 
the confining potential and to vary the electron occupation by adjusting external parametersj^l thus allowing direct 
investigation of fundamental properties of charge correlation. 

At the same time, however, these results imply that many body effects cannot be neglected in the theoretical descrip- 
tion of QDs, and may be very relevant for predictions on their possible applications, for example in semiconductor- 
based few-electron devicftSp. JjS[hile capacitance and tunneling spectroscopy experiments are becoming available for 
different dot structures Jj'LTqIj'EI it is still an open and difficult problem to devise-jtheoretical approaches allowing to 
include many body effects in calculations for realistic systems. Exact calculationaJ are necessarily limited to very few 
electrons; state-of-the-art QD structures instead involve several electrons (of the order of ten or hundred), therefore 
in general out of reach for these methods. 

In this letter we propose an accurate and effective theoretical scheme that includes electron-electron interaction 
and can be applied to systems up to a large number of electrons. This approach can treat arbitrary three dimensional 
(3D) confinement potentials, and hence deal with realistic QD structures also in the presence of external magnetic 
field. We apply this scheme to structures that were investigated in a recent experiment. The results are in good 
quantitative agreement with the observed spectra,B and provide a transparent interpretation of their physical origin. 

A key quantity that characterizes transport into a QD is the addition energy, i.e. the energy A{N) required in 
order to place an extra electron into a dot that is initially occupied by — 1 particles. Such quantity, analogous to 
electron affinity in atomic physics, can be measured experimentally as a function of A'^. A systematic experiment on 
disk-shaped dotsO has shown that the voltage increment AA between successive single electron tunneling processes 
- i.e. between two successive maxima in the-conductance - peaks at "magic" values of A^ (A^ = 2,6,12), consistent 
with the filling of complete shells calculatedQ for a dot of the same symmetry. Moreover, the experimental analysis 
showed that unusually large values of A A also occur for values of A^ cofresppnding to half-shell filling (A^ = 4, 9). The 
origin of these features, reminiscent of Hund's rule in atomic physicsJj'BllJ is intimately related to electron-electron 
interaction; therefore their quantitative evaluation is a challenge for any theoretical analysis of addition spectra in 
QDs. 

Our description is based on an accurate treatment of the single particle Hamiltonian and on the inclusion of Coulomb 
correlation according to first order perturbation theory: assuming the many particle ground state to be described by 
a single Slater determinant, the total energy of the full interacting Hamiltonian is 

E{N) = '^Eaiflaa) + ^ X! P^Pi^P-'^) + ~ -^Q/s) ("-/3(t)] (^Vct), (1) 

QO" Q/3cr 

where {fiaa) is the orbital occupation number, a denoting the set of radial and angular quantum numbers (n, m) and 
cr the spin value. According to (||), the evaluation of the energies for the various many electron states in the dot 
requires only the knowledge of the Coulomb and exchange matrix elements 
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Here k is the dielectric constant and (pa are single particle real-space eigenstates of the dot in the presence of the 
applied magnetic field. 

The important feature of the present scheme is that Uafi, JafS are calculated directly, contrary to other common 
approaches where Coulomb and exchange integrals are taken as input parameters .EJ Thus, many body effects are 
taken into account through and quantum-confinement effects through the single particle ingredients, energies 
Ea and wavefunctions (pa] the latter enter only in the calculation of Coulomb and exchange matrix elements. This 
has the advantage that realistic confining potentials (e.g. extracted from high-resolution TEM micrographs) and an 
external magnetic field can be included quite easily by chosing the appropriate single particle Hamiltonian, as is now 
commonly done for the calculation of optical properties .113 

Note that previaus, approaches have modeled the QD confinement in terms of a purely two dimensional (2D) 
coordinate spaceEoE^ As a matter of fact, in most state-of-the-art QD structures the confinement is far from being 
2D: for example, in the gated dots of Ref. |^ the typical confinement width in the growth direction z (of the order 
of lOnm) is comparable with the typical extension of the carrier ground state in the lateral-confinement (xy) plane. 
Therefore, in our calculations we fully retain the 3D nature of the problem and mimic the dots of Ref. |^ by assuming 
a finite-barrier quantum- well potential along z and a parabolic potential in (xy). Note that the assumption of 
parabolic potentials reflects the experimental indications due to the specific sample structure and is not required by 
our model. Within the usual envelope-function approximation, the only input parameters are the electron effective 
masses (to* — O.OGStoq in the dot and to* = 0.079too in the barrier), the dielectric constant in the dot (k — f2.98), 
and the quantum-well height and width {Vq = 200 meV and L = 12 nm). 

Our results for the addition energy variations, AA{N) — A{N -f 1) — A{N), are displayed in Fig. |l| as a function of 
the electron number N for different dots, characterized by different in-plane confinement energies hujQ. Here, A(N) 
is obtained as E{N) — E{N — 1), where E{N) is the ground state energy in Eq. (Q). As we can see, AA{N) exhibits 
peaks corresponding both to complete and half shell filling, and is in excellent agreement with experiments in Ref. |^. 

The origin of this result is in the delicate interplay between single particle contributions and electron-electron 
repulsion, which according to Eq. (|l]) defines the ground state configuration: the single particle term favours complete 
shell filling, while the repulsion among parallel-spin electrons, smaller than the repulsion among opposite-spin ones, 
makes the configurations with maximum total spin energetically favoured. Hund's rule is therefore already contained 
in Eq. (|l|) and clearly explains the physical origin of the half shell-filling structures. In fact, adding an electron to a 
half-filled shell forces the double occupancy of a level: consequently, AA is raised by the dominant Coulomb repulsion 
Uaa between opposite-spin electrons on the same level. 

Ground state configurations and filling rules change when a magnetic field is applied. It affects both single particle 
energies and U- and J-integrals through the induced changes in the wavefunction localization. For sufficiently large 
values, the magnetic field energetically favours configurations with higher total angular quantum number to. This is 
the physical origin of the wiggles in the A{N) vs. B plot shown in Fig. |^ and observed in the experiments of Ref. ||. 

Figure ^ shows the Coulomb and exchange integrals vs. B for the first states, obtained for Tiojq = 7.5 meV. For 
comparison, we also show the corresponding quantities calculated within the simplified 2D model. We clearly see 
that {/-integrals describing the interaction between opposite-spin electrons are smaller by a few meV in the case of 
3D confinement, while the differences in the interaction between parallel-spin electrons are much smaller. This affects 
dramatically the energy balance which determines ground state configurations, thus clearly showing the inadequacy 
of a pure 2D description of state-of-the-art QD structures. 

The proposed approach shares with the Hartree-Fock (HF) method the ansatz for total energy, expressed as the 
average of the exact Hamiltonian over a single Slater determinant; the variational prescription - allowing the construc- 
tion of optimal single particle orbitals through self-consistency - is not contained in our approach. We notice however 
that in the present case the results of perturbation theory are already much closer to the exact results than in previous 
2D calculations :I13 this is due to the fact that Coulomb integrals are artificially enhanced in 2D with respect to the 
realistic 3D description. Hence the role of self-consistency and the approximations of HF theory are much less relevant 
in 3D. Finally, we point out that the present formalism can be reduced to a Hubbard model description: expression 
(|l|) is in fact formally equivalent to the exact solution of a single-site Hubbard Hamiltonian. This is important in 
view of a direct extension of this approach to coupled dots. 

In summary, we have proposed an effective theoretical approach for the analysis of addition spectra of quantum dots 
by combining a careful 3D description of electron-electron interaction with the realistic treatment of single particle 
properties. By construction, the proposed scheme is ideally suited for applications to dots with confinement potentials 
of arbitrary shapes and symmetries. A full 3D description of Coulomb interaction is found to be crucial in determining 
single electron conductance properties of realistic nanostructures. 
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FIG. 1. Calculated addition energy increment Ayl(iV) vs. electron number N for different QD structures characterized by a 
parabolic potential in the xy plane (confining energy TvjOo) and by a finite-barrier quantum well potential along the z direction. 

FIG. 2. Addition energy A{N) vs. magnetic field calculated for a realistic (3D) QD structure with confinement energy 
Tiuio = 7.5 meV and for different electron numbers. Ground state configurations, varying with B, are also shown. 

FIG. 3. Comparison between Coulomb and exchange integrals Uap and J^p calculated assuming a simple two dimensional 
parabolic model (2D) and the corresponding three dimensional model (3D), a and (3 denote the sets of radial and angular 
quantum numbers (n, m) for the two single particle states involved in the interaction. 
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Fig. 2 M. Rontani, F. Rossi, F. IVIanghi, and E. IVIolinari 



Fig. 3 M. Rontani, F. Rossi, F. IVIanglii, and E. IVIolinari 
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